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Abstract 

Contrary to widespread perception, there is ever since 1994 a uni- 
fied, general type independent theory for the existence of solutions 
for very large classes of nonlinear systems of PDEs. This solution 
method is based on the Dedekind order completion of suitable spaces 
of piece-wise smooth functions on the Euclidean domains of definition 
of the respective PDEs. The method can also deal with associated 
initial and/or boundary value problems. The solutions obtained can 
be assimilated with usual measurable functions or even with Hausdorff 
continuous functions on the respective Euclidean domains. 
It is important to note that the use of the order completion method 
does not require any monotonicity condition on the nonlinear systems 
of PDEs involved. 

One of the major advantages of the order completion method is that 
it eliminates the algebra based dichotomy "linear versus nonlinear" 
PDEs, treating both cases with equal ease. Furthermore, the order 
completion method does not introduce the dichotomy "monotonous 
versus non-monotonous" PDEs. 

None of the known functional analytic methods can exhibit such a 
performance, since in addition to topology, such methods are signifi- 
cantly based on algebra. 



1. A Sample of customary perception 

The 2004 edition of the Springer Universitext book "Lectures on PDEs" 
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by V I Arnold, starts on page 1 with the statement : 

"In contrast to ordinary differential equations, there is no 
unified theory of partial differential equations. Some equa- 
tions have their own theories, while others have no theory 
at all. The reason for this complexity is a more compli- 
cated geometry ..." (italics added) 

The 1998 edition of the book "Partial Differential Equations" by L C 
Evans, starts his Examples on page 3 with the statement : 

"There is no general theory known concerning the solv- 
ability of all partial differential equations. Such a theory 
is extremely unlikely to exist, given the rich variety of phys- 
ical, geometric, and probabilistic phenomena which can be 
modeled by PDE. Instead, research focuses on various par- 
ticular partial differential equations ..." (italics added) 

And yet, in 1994, in Oberguggenberger & Rosinger, MR 95k:35002, 
precisely such a general theory of existence of solutions for very large 
classes of nonlinear PDEs was published. For latest developments, see 
Anguelov & Rosinger. 

In the sequel, we present the main ideas and motivations which under- 
lie the order completion method. The detailed mathematical develop- 
ments can be found in the references, among them, in Oberguggen- 
berger & Rosinger and Anguelov & Rosinger. 

It is on occasion worth recalling that we all do mathematics based on 
certain underlying ideas and motivations. What happens is that we 
may hold to them for longer, and do so long enough, so that many of 
them may become rather automatic. And once that happens, we do 
no longer - and in fact, can no longer - review them, and do so at least 
now and then. This is, then, how perceptions are established, and we 
end up being subjected to them. 

Here an attempt is made to go beyond such perceptions in the realms 
of solving PDEs. And since perceptions are inevitably formulated in a 
"meta-language" - in this case "meta" with respect to the usual formal 
mathematical texts - much of what follows has to go along with that. 
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2. The class of nonlinear systems of PDEs solved 

In Oberguggenberger & Rosinger it was show how to obtain solutions 
U for all systems of nonlinear PDEs with associated inital and/or 
boundary value problems, where the equations are of the form 

(2.1) F(x,U(x), . . . ,D?U(x), ...) = f(x), x E n C R", \p\ < m 

Here F is any function jointly continuous in all its arguments, the 
right hand term / can belong to a class of discontinuous functions, 
the order m G N is given arbitrary, while the domain Q can be any 
bounded or unbounded open set in R n . 

In fact, even the functions F defining the nonlinear partial differential 
operators in the left hand terms of (2.1) can have certain types of 
discontinuities. 

The solutions U obtained in Oberguggenberger & Rosinger can be as- 
similated with usual measurable functions on the respective domains 

a 

Recently, in Anguelov & Rosinger this general regularity result was 
further improved as it was shown that the solutions U can be associ- 
ated with Hausdorff continuous functions on the same domains Q. 

Here it is important to note that Hausdorff continuous functions are 
not much unlike usual real valued continuous functions. Indeed, on 
suitable dense subsets of their domains of definition, Hausdorff con- 
tinuous functions have as values real numbers, and are completely 
determined by such values. On the rest of their domains of definition, 
Hausdorff continuous functions can have values given by bounded or 
unbounded closed intervals of real numbers. Also, every real valued 
function which is continuous in the usual sense will be Hausdorff con- 
tinuous as well. 

One of the major advantages of the order completion method is that 
it eliminates the dichotomy between linear, and on the other hand, 
nonlinear PDEs, treating both cases with equal ease. Indeed, the di- 
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chotomy between linear and nonlinear follows from the vector space 
structure of the spaces of functions on which the partial differential op- 
erators act. In this way, this dichotomy is of an algebraic nature. On 
the other hand, partial orders are more basic mathematical structures, 
and as such, they do not, and simply cannot, differentiate between lin- 
ear and nonlinear. 

Clearly, functional analytic methods, which rely not only on topolog- 
ical but also algebraic structures cannot exhibit such a performance, 
since are bound to discriminate between linear and nonlinear, see de- 
tails in section 9. 



3. A short history of difficulties in solving linear and 
nonlinear PDEs 

The first general, that is, type independent existence result for solu- 
tions of rather arbitrary nonlinear systems of PDEs was obtained in 
1874, when upon the suggestion of K Weierstrass, Sophia Kovalevskaia 
gave a rigorous proof for an earlier theorem of Cauchy, published in 
the 1821, in his Course d' Analyse. This result although completely 
general as far as the type independent nonlinearities involved ar con- 
cerned, assumes however, that in the systems of PDEs of the form (2.1) 
both F and / are analytic. In addition one also assumes initial value 
problems on non-characteristic analytic hypersurfaces, while boundary 
value problems are not treated by the respective Cauchy-Kovalevskaia 
theorem. 

However, in such a highly particular situation concerning the regu- 
larity of the PDEs and the data involved, the solutions obtained are 
proved to exist always, and also to be unique and analytic. 
The problem with that classical existence, uniqueness and regularity 
result is that, typically for nonlinear PDEs, such analytic solutions 
do not - and in general, cannot - exist globally on the whole of the 
domain of the respective PDEs, but only in certain neghbourhoods of 
the analytic hypersurfaces on which the initial values are given. This 
is, therefore, not due to the specific method of proof of Kovalevskaia. 

Here however it is important to note that the failure of the existence 
of global analytic solutions is but a part of a far more general phe- 
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nomenon, since even linear, let alone nonlinear PDEs may fail to have 
smooth, or even merely classical solutions, even in the case of solutions 
of major applicative interest. After all, even linear constant coefficient 
PDEs have nonclassical solutions of particular interest, such as those 
given by Green functions. In the nonlinear case, difficulties start with 
the simple ODE given by U t = U 2 which does not have global classi- 
cal solutions either, except for the trivial solution U — 0. As for the 
nonlinear shock wave equation U t + UU x = 0, its nonclassical solutions 
are precisely those which model shocks. 

What may be interesting, and also worthwhile to note with respect to 
the mentioned Cauchy-Kovalevskaia theorem, are the following three 
facts : 

• The rigorous proof by Kovalevskaia of that theorem on solutions 
of general nonlinear systems of PDEs predates by about two 
decades the corresponding general theorem on solving systems 
of nonlinear ODEs defined by continuous expressions. Indeed, 
the existence of solutions for such ODEs was given by Charles 
Emile Picard in his 1894 Comptes Rendu Acad. Sci. Paris paper, 
where the associated Cauchy problem was solved by the method 
of successive approximations. 

• The only so called "hard" mathematics used in the proof of the 
Cauchy-Kovalevskaia theorem is the formula for the summation 
of a convergent geometric progression, the rest of the proof be- 
ing but a succession of rather elementary, even if quite involved, 
estimates of terms in power series. In this way, the proof of the 
Cauchy-Kovalevskaia theorem does not involve methods of func- 
tional analysis, and certainly it could not involve such methods 
at the time in the 1880s when it was given. On the other hand, 
the proof of the corresponding general existence result for solu- 
tions of nonlinear systems of ODEs does involve a fixed point 
argument in suitable spaces of functions which are complete in 
their respective topologies. 

• The result in the Cauchy-Kovalevskaia theorem - when consid- 
ered on its own original terms of type independent nonlinear 
generality - could not so far be improved in those very terms, 
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regardless of all the advances in functional analysis of the last 
more than a century. The only such improvement of the clas- 
sical result in the Cauchy-Kovalevskaia theorem was obtained 
in 1985, without however using functional analytic methods, see 
section 4. Indeed, when it comes to type independent nonlin- 
ear generality, the functional analytic methods used in solving 
PDEs could only bring about improvements - and often quite 
dramatic ones - only in a variety of far more particular cases 
than the type independent nonlinear generality dealt with in 
the Cauchy-Kovalevskaia theorem. In this way, in spite of more 
than one century of functional analysis, the classical Cauchy- 
Kovalevskaia theorem still remains a maximal result, except for 
its extension mentioned in section 4, which does not use func- 
tional analysis. 

In the early 1950s, soon after the introduction of the linear theory 
of distributions by L Schwartz, it was proved independently by Mal- 
grange and Ehrepreis that in case of a single PDE of the form (2.1), if 
the left hand term F is linear and with constant coefficients, while / is 
the Dirac delta distribution, then (2.1) always has a global, so called, 
fundamental solution given by a suitable Schwartz distribution. 

This rather general linear result appeared to suggest that a similar 
result could be obtained in the more general case when F in (2.1) 
is linear and with smooth coefficients. L Schwartz himself is known 
to have conjectured such a generalization, and furthermore, as it ap- 
pears, he suggested it at the time to Francois Treves as a subject for 
his doctoral thesis. 

However, in 1957, Hans Lewy showed that the rather simple linear 
first order PDE in three space variables and with first degree polyno- 
mial coefficients 

(3.1) (D x + iD y -2(x + y)D z )U(x,y,z) = f(x,y,z), 

(x,y,z) G R 3 

does not have any Schwartz distribution solutions in any neighbour- 
hood of any point in R 3 , for a large class of smooth right terms /. In 
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1967, Shapiro gave a similar example of a smooth linear PDE which 
does not have solutions in Sato's hyperfunctions. 

In the early 1960s, L Hormander gave certain necessary conditions for 
the solvability in distributions of arbitrary linear smooth coefficient 
PDEs, see Rosinger [7, pp. 37-39], [8, pp. 212-214]. 

Regarding the perception that nontrivial general, type independent re- 
sults are just about impossible to obtain related to PDEs, it is worth 
noting that the Malgrange-Ehrenpreis result on fundamental solutions 
is precisely such a nontrivial general and type independent existence 
result within the range of all linear and constant coefficient PDEs. 
The necessary condition for the existence of distributions solutions 
given by Hormander is also a nontrivial general and type independent 
result, this time within the much larger class of all linear smooth co- 
efficient PDEs. 



4. Nonlinear algebraic theory of generalized functions 

This nonlinear theory - see 46F30 in the AMS Subject Classifica- 
tion 2000 at www.ams.org/index/msc/46Fxx.html - was started in 
the 1960s, Rosinger [1-17], and is based on the construction of all pos- 
sible differential algebras of generalized functions which contain the 
Schwartz distributions. That theory has managed to come quite near 
to solving the Lewy impossibility. Yet it did not solve it completely, 
although it obtained generalized function solutions for large classes of 
linear and nonlinear PDEs. As an example, back in 1985, it obtained 
the first global existence result for the general nonlinear PDEs in the 
classical Cauchy-Kovalevskaia theorem. And the respective global so- 
lutions are analytic on the whole of the domain of the PDEs, except 
for certain closed and nowhere dense subsets, which can be chosen to 
have zero Lebesgue measure, see Rosinger [7, pp. 259-266], [8, pp. 
101-122], [9]. 

5. The order completion method 
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Surprisingly, the order completion method in solving general nonlinear 
systems of PDEs of the form (2.1) is based on certain very simple, even 
if less than usual, approximation properties, see Oberguggenberger & 
Rosinger [pp. 12-20]. To give here an idea about the ways the order 
completion method works, we mention some of these approximations 
here in the case of one single nonlinear PDE of the form (2.1). 

Let us denote by T(x, D) the left term in (2.1), then we have the basic 
approximation property : 

Lemma 5.1 

V x G Q, e > : 

3 5 > 0, P polynomial in x G R n : 



Consequently, we obtain : 

Proposition 5.1 

V e > : 

3 T t C Q closed, nowhere dense in Q, U e G C°°(f2) : 

f-e < T(x,D)P < f on Q \ T e 

Furthermore, one can also assume that the Lebesgue measure of T e is 
zero, namely 

mes (r e ) = 0. 



Let us now note that, see Anguelov, or Anguelov & Rosinger 





□ 



□ 



C°(fi) C H(fi) 
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and the set H(Q) of Hausdorff continuous functions on Q is Dedekind 
order complete. 

Consequently, we obtain the following basic result on the existence 
and regularity of solutions for nonlinear PDEs of the form (2.1) : 

Theorem 5.1 

T(x, D)U(x) = f(x), igU 

has solutions U which can be assimilated with Hausdorff continuous 
functions, for a class of discontinuous functions / on Q, class which 
contains the continuous functions on Q. 

□ 

We give here some more details related to Theorem 5.1 above. In view 
of Proposition 5.1, we shall be interested in spaces of piecewise smooth 
functions given by 

3 TcO closed, nowhere dense : 
*) u : Q\T -> R 

**) ueC l (n\T) 

where / G N. It is easy to see that we have the inclusions 

(5.2) T(x,D) C'n d (Vt) C C° d (Q) c H(fl) 

In this way, we obtain the following more precise formulation of the 
result in Theorem 5.1 on the existence and regularity of solutions : 

Theorem 5.1* 

(5.3) T(x,D)*(CZ(Q))# = (C&(n))# C H(fl) 

□ 



(s.i) cm = 



< u 
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Here (C™ d (fl))1p and (C° d (fi))# are Dedekind order completions of 
C™ d (Vt) and C° d (fi), respectively, when these latter two spaces are 
considered with suitable partial orders. The respective partial order 
on C™ d (£l) may depend on the nonlinear partial differential opera- 
tor T(x,D) in (5.2), while the partial order on C T ° td (f2) is the natural 
point-wise one at the points where two functions compared are both 
continuous. 

The operator T(x, D)# is a natural extension of the nonlinear partial 
differential operator T(x, D) in (5.2) to the mentioned Dedekind order 
completions. 

The meaning of (5.3) is twofold : 

• for every right hand term / G (C* d (fi))* in (2.1), there ex- 
ists a solution U G (C^(f2))*, and the set (C° d (f2)) # contains 
many discontinuous functions beyond those piecewise discontin- 
uous ones, see Oberguggenberger & Rosinger, 

• the solutions U can be assimilated with Hausdorff continuous 
functions on f2, see Anguelov & Rosinger. 

6. Comparison with methods in Functional Analysis 

The order completion method is a powerful alternative to the usual 
functional analytic ones, when solving linear or nonlinear PDEs. De- 
tails in this regard are presented in Oberguggenberger & Rosinger 
[chap. 12]. Certainly, the order completion method is not meant to 
replace the functional analytic ones, the latter being useful in obtain- 
ing stronger results in a large variety of particular PDEs. 
Here, we shall only mention the following. Functional analytic meth- 
ods in solving PDEs are based on the topological completion of uniform 
spaces, such a normed or locally convex vector spaces of suitably cho- 
sen functions. In this respect, the comparative advantages of the order 
completion method can shortly be formulated as follows : 

• unlike the functional analytic methods, which are geared more 
naturally to the solution of linear PDEs, the order completion 
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method performs equally well in the case of both linear and 
nonlinear PDEs, see section 9 below, 

• unlike the functional analytic methods, which face considerable 
difficulties when dealing with initial, and especially, boundary 
value problems, the order completion method performs without 
significant additional troubles in such situations, 

• the order completion method gives solutions which can be as- 
similated with usual measurable, or even Hausdorff continuous 
functions, and thus the solutions obtained are not merely distri- 
butions, generalized functions or hyperfuntions. 

As an illustration of the comparative situation in these two methods 
let us consider on a bounded Euclidean domain f2, which has a smooth 
boundary dfl, the following well known linear boundary value problem 

A U(x) = f(x), xen 

(6.1) 

U = on 90 

As is well known, for every given / e C°°(Q), this problem has a 
unique solution U in the space 

(6.2) X = { veC°°(n) | v = on dVl } 
It follows that the mapping 

(6.3) X3v i — > \\Av\\ L 2 {n) 

defines a norm on the vector space X. Now let 

(6.4) Y = C°°(fi) 

be endowed with the topology induced by L 2 (Q). Then in view of 
(6.1) - (6.4), it follows that the mapping 

(6.5) A : X -> Y 
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is a uniform continuous linear bijection. Therefore, it can be extended 
in a unique manner to an isomorphism of Banach spaces 

(6.6) A:X -> Y = L 2 {Q) 

In this way one has the classical existence and uniqueness result 
V /Gl 2 (0) : 

(6.7) 3! [/Gl : 

AU = f 

The power and simplicity - based on linearity and topological comple- 
tion of uniform spaces - of the above classical existence and unique- 
ness result is obvious. This power is illustrated by the fact that the 
set Y = L 2 (Q) in which the right hand terms / in (6.1) can now be 
chosen is much larger than the original Y = C°°(Q). Furthermore, 
the existence and uniqueness result in (6.7) does not need the a priori 
knowledge of the structure of the elements U E X, that is, of the 
respective generalized solutions. This structure which gives the regu- 
larity properties of such solutions can be obtained by a further detailed 
study of the respective differential operators defining the PDEs under 
consideration, in this case, the Laplacian A. And in the above specific 
instance we obtain 

(6.8) X = H 2 (n)f]H^(Q) 

As seen above, typically for the functional analytic methods, the gen- 
eralized solutions are obtained in topological completions of vector 
spaces of usual functions. And such completions, like for instance the 
various Sobolev spaces, are defined by certain linear partial differential 
operators which may happen to depend on the PDEs under consider- 
ation. 

In the above example, for instance, the topology on the space X ob- 
viously depends on the specific PDE in (6.1). Thus the topological 
completion X in which the generalized solutions U are found accord- 
ing to (6.7), does again depend on the respective PDE. 



12 



On the other hand, with the method of order completion we are no 
longer looking for generalized solutions, and instead, a type indepen- 
dent and universal or blanket regularity property is attained, since the 
solutions obtained can always be assimilated with usual measurable 
functions, or even with Hausdorff continuous functions. Similar to the 
functional analytic methods, however, the order completion method 
obtains the solutions in spaces which may again be related to the spe- 
cific nonlinear partial differential operators T(x, D) in the equations 
of form (2.1). 

7. Solving equations by extending their domains 
of definition : the three classical methods 

The ancient case of solving an equation, which shocked Pythagoras 
two and a half millennia ago, is given by 

(7.1) x 2 = 2 

This is of the general form 

(7.2) E(x) = c 

where we are given a mapping 

(7.3) E : X -> Y 

together with a specific c G Y, and then we want to find a solution 
x E X so that (7.2) holds. 

What shocked Pythagoras was that (7.1) could not be solved if one re- 
stricted oneself to X = Q in (7.3). And it took no less than about two 
millennia or more, till we could rigorously extend X = Q to X = R, 
and thus obtain a well defined solution x = ±y/2 of (7.1). 

In this way, ever since, we have the following model lesson in solving 
equations : 
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if one cannot solve (7.2) within the framework of (7.3), then one 
can try to solve it in the extended framework 



(7.4) E : X -> Y 

where X C X, Y C F, and J? is such that we have the commutative 
diagram 

x k y 



(7.5) 



c 



c 



y 



Here however, we face the following problems : 



• how to choose or construct X, and then how to interpret the new, 
or so called generalized solutions x e X\X, which two questions 
altogether constitute but the celebrated regularity problem, 

• how to do the same for Y, which nevertheless need not always 
be done, since we can often stay with Y in (7.4) and only have 
to extend X to X, 

• how to define the extension E, which often, and typically in the 
nonlinear case, is not a trivial problem. 

Fur further detail, we can now recall that with the equation (7.1) we 
had to 

(7.6) go from Q to Q = R 
On the other hand, with the equation 
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(7.7) x 2 + 1 



= 



we had to 

(7.8) go from R to C 

However, there is a vast difference between (7.6) and (7.8), respec- 
tively, between solving (7.1) and (7.7). Indeed, we solve (7.7) through 
the extension (7.8) which is a mere algebraic adjoining of an element, 
in this case, of i = to R. 

On the other hand, when solving (7.1), the extension (7.6) can be seen 
as at least three different, even if in this particular case equivalent, 
constructions, namely, through : 

• topology 

• algebra 

• order. 

And to be more precise, we have : 

• The Cauchy-Bolzano method is ring theoretic plus topological, 
and it is applied to Q, as it obtains R according to the quotient 
construction in algebras 

(7.9) R = A/I 

where A C Q N is the algebra of Cauchy sequences of rational 
numbers, while I is its ideal of sequences convergent to zero. 

• The method of Dedekind is based on the order completion of Q. 

However, the Cauchy-Bolzano method can be generalized in two di- 
rections : 

• In the topological generalization the algebraic part can be omit- 
ted, and instead, one only uses the topological completion of uni- 
form spaces, here of the usual metric space on Q. 
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• In the algebraic generalization it is possible to extract the ab- 
stract essence of (7.9), and simply start with a suitable algebra 
A, and an ideal X in it. Such a construction can indeed be rather 
abstract, since it need not involve any topology on Q or R, as it 
happens for instance, when constructing the nonstandard reals 
*R, namely 

*R = A/1 

Here one takes A = R N , that is, the algebra of all sequences 
of real numbers, while the ideal X is defined by any given free 
ultrafilter on N. 

A rather general version of such an abstract approach, which 
however makes a certain limited use of topology, has been intro- 
duced and extensively used in the nonlinear algebraic theory of 
generalized functions under the AMS classification index 46F30, 
as mentioned in section 4 above. 

What is done in the method in 46F30 is to generalize the Cauchy- 
Bolzano method by retaining its ring theoretic algebraic aspect, while 
the topological one is weakened to the certain extent of being confined 
to the topologies of Euclidean spaces only. 

What is done in the method introduced in Oberguggenberger & Rosinger, 
and further developed in Anguelov & Rosinger, is the extension of the 
classical Dedekind order completion method, used in the construction 
of R from Q, to suitable spaces of piece-wise smooth functions. 

An important fact to note is that both the topological and order com- 
pletion methods give us the property that 

• Q is dense in R 

in the respective sense of topology or order. In this way, the elements 
of the extension of Q, that is, the elements of R, are in the corre- 
sponding sense arbitrarily near to the elements of the extended space 
Q. Thus the elements in the extension can arbitrarily be approximated 
by elements of the extended space, be it in the sense of topology, or 
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respectively, order. 

Furthermore, both through the methods of topology and order, one 
obtains R in a unique manner, up to a respective isomorphism. 

In this way both the topological and order completion methods have 
the double advantage that 

• the elements of the extension are not too strange conceptually, 
and furthermore 

• the elements of the extension are near to elements of the ex- 
tended space, within arbitrarily small error. 

This density property remains also in the general Dedekind order com- 
pletion method used in Oberguggenberger & Rosinger and Anguelov 
& Rosinger. Indeed, in (5.3) we have that C^(Q) and C° d (fi) are 
order dense in {C™ d {Vt))* and (C° d (f2)) # , respectively. 

Connected with the general extension method in (7.3) - (7.5) one can 
note that, on occasion, the following convenient situation may occur : 
the extended mapping 

(7.10) E:X — ► Y 

may turn out to be an isomorphism of the respective algebraic, topo- 
logical or order structures used on X and Y, when constructing the 
corresponding extensions X and Y. In such and when one 

has a better understanding of the structure of the elements in Y, one 
can obtain in addition a regularity type result concerning the so called 
generalized solutions x E X \ X of the equations (7.2), since such 
generalized solutions can be assimilated - through the isomorphism E 
- with the corresponding elements E(x) G Y. 

A classical example of such an isomorphism (7.10) happens, for in- 
stance, in (6.5), (6.6), when the boundary value problem (6.1) is solved 
by using well known functional analytic methods. 
In that specific instance, however, the suitable further use of func- 
tional analytic methods can lead to the additional regularity property 
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of generalized solutions in X, as given in (6.8). Nevertheless, the Ba- 
nach space isomorphism (6.6) - which in that case is but the particular 
form taken by (7.10) - is in itself already a first regularity result about 
the structure of the elements of X. 

The above convenient situation of an isomorphism of type (7.10) can 
appear as well when using the order completion method in solving non- 
linear system of PDEs. This is the reason why the solutions obtained 
in Oberguggenberger & Rosinger could be assimilated with usual mea- 
surable functions, while in Anguelov & Rosinger, they can be assimi- 
lated with the much more regular Hausdorff continuous functions. 
More specifically, in (5.3), the extended mappings T(x, -D) # prove to 
be order isomorphisms between the spaces C^(f2)* and C° d (f2) # . 
This is then, in essence, the reason why the solutions of nonlinear 
systems of PDEs of the form (2.1) could earlier be assimilated with 
usual measurable functions, and can now be assimilated with Haus- 
dorff continuous functions. 



8. The need for extensions in the case of solving PDEs. 

Let us now associate with each nonlinear PDE in (2.1) the correspond- 
ing nonlinear partial differential operator defined by the left hand side, 
namely 

(8.1) T(x,D)U(x) = F(x,U(x), . . . ,D>U(x), . . . ), xeVl 

Two facts about the nonlinear PDEs in (2.1) and the corresponding 
nonlinear partial differential operators T(x,D) in (8.1) are important 
and immediate 

• The operators T(x, D) can naturally be seen as acting in the 
classical context, namely 

(8.2) T(x, D) : C m (n) 3 U i — ► T(x, D)U G C°(ft) 

while, unfortunately on the other hand 

• The mappings in this natural classical context (8.2) are typi- 
cally not surjective. In other words, linear or nonlinear PDEs in 
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(2.1) typically cannotbe expected to have classical solutions U G 
C m (Q), for arbitrary continuous right hand terms / G C°(fl). 
Furthermore, it can often happen that nonclassical solutions 
have a major applicative interest, thus they have to be sought 
out beyond the classical framework in (8.2). 

This is, therefore, how we are led to the necessity to consider general- 
ized solutions U for PDEs of type (2.1), that is, solutions U £ C m (f2), 
which therefore are no longer classical. This means that the natural 
classical mappings (8.2) must in certain suitable ways be extended to 
commutative diagrams 



C m (tt) 



T(x,D) 



C°(fi) 



i.3) 



C 



c 



X 



y 



T 



which are expected to have certain kind of surjectivity type properties, 
such as for instance 



(8.4) c°(n) C T(X) 

We conclude with a few comments : 



• Traditionally, ever since Hilbert and Sobolev, starting before 
WW II, functional analysis has been used in solving PDEs, 
and suitable uniform topologies are defined on the domains and 
ranges of the corresponding partial differential operators T(x, D). 
Thus these operators obtain certain continuity properties. Then 
the extensions X and y in (8.3) are defined as the completions in 
these uniform topologies of the domains and ranges of T(x, D), 
respectively. Finally, the continuity properties of T(x, D) may 
allow the construction of suitable extensions T which would give 
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the commutative diagrams (8.3), and also satisfy some version of 
the surjectivity property (8.4), see for details Oberguggenberger 
& Rosinger [chap. 12, pp. 237-262]. 

• Since the 1960s, the algebraic method in 46F30 can alternatively 
be used especially in the case of nonlinear partial differential 
operators T(x,D). In this respect, large classes of differential 
algebras of generalized functions containing the Schwartz distri- 
butions were constructed as the sought after extensions X and 
y in (8.3). 

The most general classes of such algebras were introduced and 
used in Rosinger [1-17], starting with the 1960s. Later, in the 
1980s, a particular class of such algebras was introduced in 
Colombeau, and it has known a certain popularity. However, 
due to the specific polynomially limiting growth conditions re- 
quired in the construction of Colombeau algebras, their use in 
the study, for instance, of Lie group symmetries of PDEs, or 
singularities in General Relativity is limited, since in both cases 
one may have to deal with transformation whose growth can be 
arbitrary. In this way, such transformations cannot be accommo- 
dated within the Colombeau algebras of generalized functions. 
On the other hand, arbitrary smooth transformations and op- 
erations can easily be dealt with in the much larger classes of 
algebras of generalized functions introduced earlier in Rosinger 
[1-17]. 

• The order completion method, introduced and developed in 1994 

in Oberguggenberger & Rosinger, and further improved in Anguelov 
& Rosinger, constructs the extensions X and y in (8.3) as the 
Dedekind order completion of spaces naturally associated with 
the partial differential operators T(x, D), and the spaces C m (Q) 
and C°(n) in (8.2). 

Related to the advantages of the order completion method in solving 
nonlinear PDEs let us mention here in short the following. 

Neither the functional analytic, nor the algebraic methods can so far 
come anywhere near to solve nonlinear PDEs of the generality of those 
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in (2.1), let alone systems of such nonlinear PDEs together with asso- 
ciated initial and/or boundary value problems. 

In fact, the functional analytic methods are still subjected to the cele- 
brated 1957 Hans Lewy impossibility which they are nowhere near to 
manage to overcome, even in the general smooth coefficient linear case. 

As far as the algebraic method is concerned, it has among others come 
quite near to the solution of the Lewy impossibility, see Colombeau, 
and for a short respective account Rosinger [7]. 

Further powerful results in solving various classes of nonlinear PDEs, 
not treated so far by the functional analytic method, can be found in 
Rosinger [5-10], Colombeau , or Oberguggenberger. 
Among such results is the global solution of arbitrary analytic nonlin- 
ear systems of PDEs, when considered with analytic non-characteristic 
Cauchy initial values, mentioned in section 4 above. The respective 
generalized solutions obtained are analytic functions, except for closed 
and nowhere dense subsets T of the domains Q of definition of the 
given PDEs. In addition, these subsets T can also be chosen to have 
zero Lebesgue measure. 

On the other hand, the order completion method introduced and de- 
veloped in Oberguggenberger & Rosinger, and further improved in 
Anguelov & Rosinger, as far as the regularity of solutions is concerned, 
can not only deliver global solutions for systems of nonlinear PDEs of 
the generality of those in (2.1), but it can also obtain a universal or 
blanket type independent regularity result for such solution, namely, 
prove that they can be assimilated with usual measurable functions, 
or even with Hausdorff continuous functions. 

Clearly, as one of the consequences of solving nonlinear systems of 
PDEs of the generality of those in (2.1), the order completion method 
in Oberguggenberger & Rosinger and in Anguelov & Rosinger is the 
only one so far which fully manages to overcome the Lewy impossibil- 
ity. Furthermore, it does so with a large nonlinear margin. 

9. Order Completion Abolishes the Dichotomy "Linear 
versus Nonlinear" 
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The dichotomy linear versus nonlinear relating to equations or opera- 
tors is in its essence an issue of algebra, and more specifically, of vector 
space structures. In this way, it is present both in the functional an- 
alytic and algebraic methods for solving PDEs. And needless to say, 
dealing with the nonlinear case proves to be incomparably more diffi- 
cult than it is with the linear one. Consequently, and unfortunately, 
the presence of this dichotomy is one of the major disadvantages of 
both the functional analytic and algebraic methods in solving PDEs, 
even if by now it is taken so much for granted that no attempt is made 
to abolish it. 

On the other hand, order structures are of a more basic type than the 
algebraic ones. 

Consequently, if instead of algebraic structures we consider order struc- 
tures on the spaces of smooth functions on which the partial differ- 
ential operators act naturally when seen in the classical context, see 
for instance (8.2), then these order structures - being more basic than 
algebra - can no longer distinguish between the linearity or nonlinear- 
ity of such partial differential operators. In this way, the traditional 
dichotomy between the linear and nonlinear is simply set aside, and 
then the only problem left is whether indeed one can solve PDEs in 
the completion of such order structures. 

Fortunately, as shown in Oberguggenberger & Rosinger, and also in 
Anguelov & Rosinger, such a solution is possible and useful. 

Here of course, one may think that it may help if the respective par- 
tial differential operators are monotonous. And then one may be con- 
cerned that all we managed to do was simply to get rid of the di- 
chotomy between linear and nonlinear, so that instead, now we have 
to face the dichotomy monotonous versus arbitrary partial differential 
operators, with the latter being quite likely again far more difficult to 
deal with, than the former. 

Such a particular approach in which the dichotomy monotonous ver- 
sus non-monotonous prevails has recently been pursued in Carl & 
Heikkila, for instance. 

This however is clearly not the way in the method in Oberguggen- 
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berger & Rosinger and in Anguelov & Rosinger Indeed, although in 
this method order structures and order completions are essentially 
used in solving PDEs, one need not require a priori any sort of mono- 
tonicity property related to the equations solved. Certainly, the gener- 
ality of nonlinear PDEs in (2.1) or of systems of such nonlinear PDEs 
illustrates the fact that the respective equations are not supposed to 
satisfy any a priori monotonicity conditions whatsoever. 

What happens is very simple in fact, and is similar with the way 
the operators T(x, D) and T in (8.3) acquire continuity type proper- 
ties when the functional analytic method is used in the construction 
of such commutative diagrams. Indeed, with the functional analytic 
method, when one starts, say, with the natural mappings (8.2), the 
uniform topologies one considers on the classical domains and ranges 
of the operators in order to obtain commutative diagrams (8.3) are 
not arbitrary, but typically are related to the respective operators, see 
section 6 above, or Oberguggenberger & Rosinger [chap. 12]. 

The same happens when order completion is used in Oberguggenberger 
& Rosinger and in Anguelov & Rosinger for the construction of com- 
mutative diagrams (8.3). More precisely, the order structures on the 
spaces of smooth functions on which the partial differential operators 
naturally act, see for instance (8.2), will typically be defined depen- 
dent to a certain extent on these operators. In this way, the respective 
operators, no matter how arbitrary within the class of those in (2.1), 
will nevertheless become monotonous, thus so much easier to deal with. 

Such a procedure obviously cannot be imitated within algebra, since 
a nonlinear operator cannot in general define a vector space structure 
in which it would become linear. 

Furthermore, such a procedure obviously goes far beyond the approach 
in Carl & Heikkila, for instance, where one starts with given natural 
order relations on both the domains and ranges of ODEs or PDEs, 
and then severely restricts oneself only to those rather small classes 
of equations whose associated operators, or rather, inverse operators, 
are monotonous in the a priori given orders. 



23 



10. The Hidden Power of Methods Based on Partial Orders 

As it happens, there is a rather widespread perception in mathematics 
that order structures are too simple, and thus powerless, especially in 
analysis, therefore, they can deliver less than algebra, which on its 
turn, can deliver less than topology. 

Accordingly, since the emergence of functional analytic methods in 
the solution of PDEs at the beginning of the 20th century, with the 
respective wealth of topologies on a variety of spaces of functions, the 
perception prevails that there simply cannot be any other more pow- 
erful methods in present day mathematics which could deal with such 
equations. 

In view of such a perception it may appear surprising to see results 
such as in Anguelov & Rosinger, and the earlier ones in Oberguggen- 
berger & Rosinger, results obtained through order structures, and 
which so far could not be approached anywhere near by functional 
analytic methods. 

Indeed, this method - based on order completion - does solve systems 
of PDEs of the nonlinear generality of those in (2.1), together with 
associated initial and/or boundary value problems, and furthermore, 
delivers for them global solutions which can be assimilated with usual 
measurable or even Hausdorff continuous functions. It is in this way 
that the order completion method is not only unprecedented, but it 
may also look rather strange in view of the mentioned perception in 
mathematics related to order structures. 

Therefore, one should indeed address the apparent secret of the power 
of the order completion method in solving such large classes of sys- 
tems of nonlinear PDEs, together with their associated initial and/or 
boundary value problems. And one can do so by questioning the men- 
tioned perception. This can perhaps best be done by the presentation 
of certain classical - even if less well known - examples which illustrate 
the power of order structures in yielding what usually are called deep 
theorems. 

In this regard a rather impressive, yet less well known fact is given 
by the 1936 "Spectral Theorem" of Freudenthal, see Luxemburg & 
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Zaanen [chap. 6]. 

Let us recall here in short some of its rather deep consequences. 
The mentioned "Spectral Theorem" is a theorem about partially or- 
dered structures, and it was proved by Freudenthal exclusively in terms 
of such structures. Yet, what is of special relevance is that by suitable 
particularizations, one can obtain from it the following three results 
which are in fields as diverse as Operator Theory, Measure Theory 
and linear PDEs : 

• the celebrated spectral representation theorem for normal oper- 
ators in Hilbert spaces, 

• the highly nontrivial Radon-Nikodym theorem in measure the- 
ory, and 

• the Poisson formula for harmonic functions in an open circle. 



11. Conclusions 

The unprecedented power of the order completion method in solv- 
ing very general systems of nonlinear PDEs and the associated initial 
and/or boundary value problems stems from two facts : 

• Partial orders are more basic mathematical structures than al- 
gebra or topology. And being more basic than algebra, partial 
orders do not distinguish between linear and nonlinear equa- 
tions, operators, and so on. Consequently, partial orders treat 
the linear and nonlinear cases in the same manner. Functional 
analytic method clearly cannot do the same. 

• When using order completion for solving PDEs, one need not 
assume any monotonicity properties of the respective equations. 

In the order completion method, the partial orders on the spaces of 
functions which are the domains of definition of the partial differen- 
tial operators considered are defined in relation to these operators. 
This is similar to the way the topologies on such domains are defined, 
when functional analytic methods are used. Further details in this 



25 



regard can be found in Oberguggenberger & Rosinger [chap. 12, 13], 
or Anguelov & Rosinger. 
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